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Abstract 

We construct a dark energy model where a scalar field non-minimally coupled to 
gravity plays the role of the dark component. We compare cosmological consequences of 
this non-minimal coupling of the scalar field and gravity in the spirit of the dark energy 
paradigm in Jordan and Einstein frames. Some important issues such as phantom 
divide line crossing, existence of the bouncing solutions and the stability of the solutions 
are compared in these two frames. We show that while a non-minimally coupled scalar 
field in the Jordan frame is a suitable dark energy component with capability to realize 
phantom divide line crossing, its conformal transformation in the Einstein frame has 
not this capability. The conformal transformation from Jordan frame to Einstein frame 
transforms the equation of state parameter of the dark energy component to its minimal 
form with a redefined scalar field and in this case it is impossible to realize a phantom 
phase with possible crossing of the phantom divide line. 
PACS: 98.80.-k, 98.80.Cq 
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1 Introduction 



Recently, evidences from supernova searches data [1,2], cosmic microwave background (CMB) 
results [3-5] and also Wilkinson Microwave Anisotropy Probe (WMAP) data [6,7], indicate 
an positively accelerating phase of the cosmological expansion today and this feature shows 
that the simple picture of universe consisting of a pressureless fluid is not enough. In this 
regard, the universe may contain some sort of additional negative-pressure dark energy. 
Analysis of the five year WMAP data [8-11] shows that there is no indication for any signifi- 
cant deviations from Gaussianity and adiabaticity of the CMB power spectrum and therefore 
suggests that the universe is spatially flat to within the limits of the observational accuracy. 
Further, the combined analysis of the three-year WMAP data with the supernova Legacy 
survey (SNLS) [8], constrains the equation of state w^e, corresponding to almost 74% contri- 
bution of dark energy in the currently accelerating universe, to be very close to that of the 
cosmological constant value. Moreover, observations appear to favor a dark energy equation 
of state, Wde < — 1 [12,13]. Therefore a viable cosmological model should admit a dynam- 
ical equation of state that might have crossed the value Wd e = —1, in the recent epoch of 
cosmological evolution. In fact, to explain positively accelerated expansion of the universe, 
there are two alternative approaches: incorporating an additional cosmological component 
to the matter part of the Einstein equations or modifying gravity ( the geometric part of 
the Einstein field equations) at cosmological scale. Multi-component dark energy with at 
least one non-canonical phantom field is a possible candidate of the first alternative. This 
viewpoint has been studied extensively in literature ( see [14] and references therein ). 

On the other hand, despite of several successes, the standard model of cosmology suffers 
from a series of problems. The most serious of these problems is the problem of initial 
singularity because the laws of physics break down at the singularity point. In order to avoid 
this lawlessness, there is a huge interest in the solutions that do not display divergencies. 
These solutions could be obtained at a classical level or by quantum modifications. Most 
of the efforts in quantum gravity is devoted to reveal the nature of the initial singularity 
and to understand the origin of matter, non-gravitational fields, and the very nature of 
the spacetime. In recent analysis done within the loop quantum cosmology, the Big Bang 
singularity is replaced by a quantum Big Bounce with finite energy density of matter. This 
scenario has strong quantum effects at the Planck scale too. Another motivation to remove 
the initial singularity is the initial value problem. A sound gravitational theory needs to 
have a well-posed Cauchy problem. Due to the fact that the gravitational field diverges at 
the singularity, we could not have a well-formulated Cauchy problem as we cannot set the 
initial values at that point. 

With these preliminaries, the purpose of the present paper is to study some currently 
important cosmological issues such as phantom divide line crossing, avoiding singularities by 
realization of the bouncing solutions and the stability of these solutions in a non-minimally 
coupled scalar field model of universe in Jordan and Einstein frames and in the spirit of dark 
energy model. We compare cosmological consequences of the non-minimal coupling between 
scalar field and gravity in the spirit of the dark energy scenario in these two frames. Some 
important issues such as phantom divide line crossing, existence of the bouncing solutions 
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and the stability of the solutions are compared in these two frames. Especially, we analyze 
the parameter space of the model numerically to show that which frame allows for stability 
of the solutions in the separate regions of the u> — u phase-plane. We show that with only 
one scalar field non-minimally coupled to gravity, crossing of the phantom divide line can 
be realized just in the Jordan frame. By transforming to Einstein's frame, we show that 
this model cannot account for crossing of the phantom divide line. We show that while a 
non-minimally coupled scalar field in the Jordan frame is a suitable dark energy component 
with capability to realize phantom divide line crossing, its conformal transformation in the 
Einstein frame has not this capability. In fact, conformal transformation from Jordan frame 
to Einstein frame transforms the equation of state parameter of dark energy component to 
its minimal form with a redefined scalar field and in this case it is impossible to realize a 
phantom phase with possible crossing of the phantom divide line. 



2 A Ricci-Coupled Scalar Field Model in the Jordan 
Frame 

For a model universe with a non-minimally coupled scalar field as matter content of the 
universe, the action in the absence of other matter sources in the Jordan frame can be 
written as follows 



= J d 4 xy/^g^a((f))R[g} - I^V^V^ - V(<f>) 



(1) 



where we have included an explicit and general non-minimal coupling of the scalar field and 
gravity as a(<f>). For simplicity, from now on we set k 2 = 8tcGn = 1- Variation of the action 
with respect to the metric gives the Einstein equations 

-^V ~~ ~^9nvR = ol 1 T fJjU . (2) 

T^ v , the energy- momentum tensor of the scalar field non-minimally coupled to gravity, is 
given by 

%v = V M 0V,0 - ^(V</>) 2 - g^V(4>) + g^Ua{<t)) - V M V,a(0), (3) 

where □ shows 4-dimensional d'Alembertian. For FRW universe with line element defined 

as 

ds 1 = -dt 2 + a 2 {t)d^ 2 , (4) 

where cE^ 2 is the line element for a manifold of constant curvature k = +1,0,-1, the 
equation of motion for scalar field (p is 

WV^ = V'-a'R[g], (5) 
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where a prime denotes derivative of any quantity with respect to (p. This equation can be 
rewritten as 

4> + 3-<P+^- = a'R[g}- (6) 

a acp 

where a dot denotes the derivative with respect to cosmic time t and Ricci scalar is given 
by 



k 



R = Q[H + 2H 2 + —y (7) 

With this non-minimally coupled scalar field as matter content of the universe, cosmological 
dynamics are described by 

or cr 3 

and 

;~ic + *>- (») 

In these equations, the effect of the non-minimal coupling of the scalar field and gravity is 
hidden in the definition of p and p. We assume that scalar field, <p, has only time dependence 
and using (3), we find 

p = a' 1 Q</> 2 + V(<f>) - 6a'H^j , (10) 

p = a' 1 Q</> 2 - V(<i>) + 2(a'4> + 2Ha'<f) + a"</> 2 )) , (11) 
where H = ^ is Hubble parameter. Now, equation (9) takes the following form 

- a = -ia- 1 ^ 2 - 2V( ( f>) + 6(a'0 + + a"^, (12) 

and dynamics of the equation of state parameter is given by 

p 2 -2V{<P)+4(a'4> + 2Ha'<j) + a f '<p 2 ) 

W = - = : " : -. (13) 

p <p 2 + 2V{<j))-l2a'H<t) 

From this equation, when = 0, we obtain p = —p. In this case p is independent of a and 
V(<p) plays the role of a cosmological constant. In the minimal case when 2 < V(<p), using 
(9) we obtain p < — | which shows an accelerated expansion which is driven by cosmological 
constant. However, cosmological constant is not a good candidate for dark energy since its 
suffers from several conceptual problems such as its unknown origin and also need to huge 
amount of fine-tuning. In non-minimal case the cosmological dynamics depends on the value 
of the non- minimal coupling. In order to solve the Friedmann equation (8), we choose the 
following scalar field potential 

V(<j>) = ^m 2 <p 2 (14) 
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With this potential, a possible solution of our basic equations, (6), (8) and (10) is as follows 
( see [15] for a similar argument) 

4> = \J~C~ocos(mi) (15) 

where Co is a parameter with the dimension of mass squared describing the oscillating 
amplitude of the fields. Here we assume a(<f>) — \{1 — £0 2 ). We note that there are two 
critical values of <fi given by (j) c = ±-K= that should be avoided to have well-defined field 

equations. Using (14) and (15) in (10), we find 

m 2 Co — 12^HComsin(mt)cos(mt) 
P= 1 - £C cos 2 {mt) ' ^ ' 

where those values of the cosmic time coordinates that lead to singular energy density are 
excluded from our considerations. Therefore, for a flat spatial geometry, Friedmann equation 
(8) can be rewritten as follows 



H 



6 cos (mt) sin (mt) C £ ± ^36 (cos (mt)) 2 (sin (mt)) 2 C 2 £ 2 + C — £C 2 (cos (mt)) 2 J m 



-1 +£C cos 2 (mt) 

(17) 

We avoid imaginary values of the Hubble parameter in which follows. 
2.1 Crossing of the phantom divide line 

Non-minimal coupling of the scalar field and gravity in the Jordan frame provides a suitable 
framework for explanation of the late-time accelerated expansion [16]. On the other hand, 
after substituting corresponding relations for <fi, H and V in equation (13), dynamics of the 
equation of state parameter for a non-minimally coupled scalar field is given by 



uj(t) = I — 8£esinmt cosmty — 36 ^— + C £ 2 cos 4 mt — £C ^£ — cos 2 mtj C 
+48£C (-^4 + f) cos4 mt + ( 2 + - 52C £ 2 ) cos 2 mt - 1 + 4£ j x 
- 12£e sin mt cos mtJ —36 ^— + C £ 2 cos 4 mt — £C ^ — cos 2 mt^j C — 1 



+72C £ 2 cos 4 mt - 72£ {-^ + ^ C cos 2 mt 



(18) 



where e = ±1. Figure 1 shows the crossing of the phantom divide line with equation of state 
parameter of this non-minimally coupled scalar field. Nevertheless, as figure 2 shows, in the 
case of £ = 0, that is a single and minimally coupled scalar field, there is no crossing of the 
phantom divide line, as has been emphasized by other literature such as [14]. We note that 
our analysis shows that for negative values of £ ( corresponding to anti-gravitation) and any 
values of e ( e = — 1 or e = +1), it is impossible to realize the phantom divide line crossing 
in this setup. 
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1 2 3 4 5 6 7 1 2 3 4 5 6 7 

t t 

| for xi > and epsilon > | | for xi > and epsilon < 0| 

Figure 1: A non- minimally coupled scalar field in Jordan frame has the capability to realize crossing 
of the phantom divide line by its equation of state parameter in a suitable subspace of the model 
parameter space. 




t 

| for xi-0 and epsilon < or epsilon > 0| 



Figure 2: Equation of state parameter of a single, minimally coupled scalar field ( with £ = ), 
cannot explain crossing of the phantom divide line. 
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2.2 Bouncing behavior of the model 

A bouncing universe with an initial contraction to a non-vanishing minimal radius and a 
subsequent expanding phase provides a possible solution to the singularity problem of the 
standard Big Bang cosmology. For a successful bounce, it can be shown that within the 
framework of the standard 4-dimensional Friedmann-Robertson- Walker (FRW) cosmology 
with Einstein gravity, the null energy condition (NEC) is violated for a period of time around 
the bouncing point. Moreover, for the universe entering into the hot Big Bang era after the 
bouncing, the EoS of the matter content w in the universe must transit from w < — 1 to 
w> -1 [17]. 

The model proposed to understand the behavior of dark energy with an EoS of w > — 1 in 
the past and w < — 1 at present, has been supported by the observational data [18]. This is a 
dynamical model of dark energy and it differs from the cosmo logical constant, Quintessence, 
Phantom, K-essence and so on in the determination of the cosmological evolution. Here we 
study the possibility of realization of the bouncing solution in a model universe dominated 
by the nonminimally coupled scaler field as matter content. We start with a detailed exam- 
ination of the necessary conditions required for a successful bounce. During the contracting 
phase, the scale factor a(t) is decreasing, i.e., a(t) < 0, and in the expanding phase we have 
a(t) > 0. At the bouncing point, a(t) = 0, and around this point a(t) > for a period of 
time. Equivalently, in the bouncing cosmology the Hubble parameter H runs across zero 
from H < to H > and H = at the bouncing point. In our model this behavior is shown 
in figure 3. A successful bounce requires that around this point 

H = -AnGp(l +w) > . (19) 

From equation (18) one can see that w < — 1 in a neighborhood of the bouncing point. After 
the bounce the universe needs to enter into the hot Big Bang era, otherwise the universe 
filled with the matter with an EoS w < — 1 will reach the big rip singularity as what happens 
to the Phantom dark energy which violates the null energy condition [19]. This requires the 
EoS of the matter to transit from w < — 1 to to > — 1. 

By solving the Friedmann equation (17), we can study the scale factor versus the cosmic 
time, t. If we integrate equation (17), we find 



ait) — aQexp 



, (mt) sin (mt) C £ ± \J 36 (cos (mt)) 2 (sin (mt)) 2 C 2 £ 2 + C a - (C 2 (cos (mt)) 2 ^j 
-1 + £C (cos (mt)) 2 



(20) 



Using equation (20), we plot the behavior of the scale factor versus the cosmic time. One 
can see from Figures 3 (left) and 4 that a nonsingular bounce happens when the Hubble 
parameter H runs across zero with a minimal and non-vanishing scale factor a. 

2.3 Stability of the model 

Now we study the stability of our model. In order to study the classical stability of our 
model, we analyze the behavior of the model in the u — u' plane where u' is the derivative 
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2 3 4 5 6 7 

t 



- for xi > and epsilon > 0| 



- for xi > and epsilon < | 



Figure 3: Variation of the Hubble parameter H relative to cosmic time t. The left hand side figure 
shows that the universe can switches alternatively between expanding and contracting phases. 




t 



Figure 4: Variation of the scale factor a relative to cosmic time t. Note that the initial conditions 
have been set so that integrand of equation (20) has non-vanishing derivative at the minimum 
point. 
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of u) with respect to the logarithm of the scale factor ( see [20-23] for a similar analysis for 
other interesting cases) 

, dui du dt lo . , 

J= = = — . (21) 

din a dt din a H 

The sound speed expresses the phase velocity of the inhomogeneous perturbations of the 
quintessence field. We define the function c a as 

c a 2 = ? (22) 

p 

If the matter is considered as a perfect fluid, this function would be the adiabatic sound 
speed of this fluid. We note that with scalar fields that do not obey perfect fluid form 
necessarily, this quantity is not actually a sound speed. In which follows, we demand that 
c a 2 > in order to avoid the big rip singularity at the end of the universe evolution. The 
conservation of the quintessence field effective energy density can be stated as 

dp quintessence ott( _i_ \ _ n (O^ 

^ ~r Oil \p quintessence ~r p 'quintessence) U 

Since the dust matter obeys the continuity equation and the Bianchi identity keeps valid, 
total energy density satisfies the continuity equation. From above equation, we have 

p de = -3Hp de (l + u de ) (24) 

Using equation of state p de = uj de p de , we obtain 

Pde = Upde + Udepde (25) 

Therefore, the function c a 2 could be rewritten as 

c ' 2= ^mr^) + ^ (26) 

In this situation, the u> — u' plane is divided into four regions defined as follows 

/: u de > -1, J > 3u(l + uS) =>• c a 2 > 

// : u de > -1, J < 3a; (1 + u) ^> c a 2 < 

/// : u de < -1, J > 3cu(l + uj) ^> c a 2 < 

IV: u de <-l, J < 3u(l + uj) c a 2 >0 (27) 

As one can see from these relations, the regions I and IV have the classical stability in our 
model. We plot the behavior of the model in the u — a/ phase plane and identify the regions 
mentioned above in figure 5. 
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-1.5 



-1.0 
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-0.5 
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Figure 5: Bounds on u/ as a function of u> in uj — uj' phase plane. The stable regions are I and IV. 

3 A Ricci-Coupled Scalar Field in the Einstein Frame 

Now we extend our study of the mentioned cosmological issues to the Einstein frame by 
adopting a conformal transformation. The action (1) in the Jordan frame can be rewritten 
as follows 



S = J d^x^g" 



2 2 ^ 



^Rcf> 2 - V{<t>) 



(28) 



where we assumed k 2 = 1 and a((f>) = |(1 — £0 2 ) and £ is a non-minimal coupling. The 

metric signature convention is chosen to be (H ) with spatially flat Robertson- Walker 

metric as follows 



ds 2 = dt 2 - a 2 (t)5 ij dx i dx j . 



(29) 



To obtain the fundamental background equations in Einstein frame, we perform the following 
conformal transformation 



g^u = Qg^, SI = 1 - i4> 2 



(30) 



Here we use a hat on a variable defined in the Einstein frame. The conformal transformation 
gives ( see for instance [16] and references therein) 



d xJ—g 



(31) 
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where by definition 



F\<t>) 



and 



v(4>) 



1 _ ^(i _ eg 

(i-£0 2 ) 2 

V(0 



(1-^2)2- 

Therefore, one may redefine the scalar field as follows 

# _ _ - e0 2 (i - 6Q 

H ~ m ~ 



(32) 



(33) 



(34) 



When we investigate the dynamics of universe in the Einstein frame, we should transform 
our coordinates system to make the metric in the Robertson- Walker form 



and we obtain 



a = V&a, di = VQdt, 



ds 2 = di 2 - a 2 (i)5 ij dx t dx j . 



(35) 



(36) 



Note that the physical quantities in Einstein frame should be defined in this coordinate 
system. Now the field equations can be written as follows 



H 2 



1 /d(j) 
2\dl 



p 

3' 



(37) 



d 2 6 n -dS dV , , 

-£ + 3H-£ + — r = 38 
dt 2 dt d<f> 

where H — ™. 

a 

We assume that scalar field has only time dependence and we find dynamics of equation 
of state as follows 

p i(#V-w) 



V* = 1 = f .(2 • ( 39 ) 

p l(d± 

di 



+ v(4>) 



This is an interesting result: while a non-minimally coupled scalar field in the Jordan frame 
is a suitable dark energy component with capability to realize phantom divide line crossing, 
its conformal transformation in the Einstein frame has not this capability. In fact, conformal 
transformation from Jordan frame to Einstein frame transforms the equation of state param- 
eter of dark energy component to its minimal form with a redefined scalar field and in this 
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Figure 6: Variation of the Hubble parameter H and scale factor relative to the cosmic time t. 
There is no possibility to realize bouncing solution in this frame. 

case it is impossible to realize a phantom phase with possible crossing of the phantom divide 
line. In a minimal coupling case with a = 0, one can obtain equation of state parameter 
cu — 1 — 2cos 2 {mi). As figure 2 shows, this minimally coupled scalar field cannot explain 
the crossing of the phantom divide line. However, we note that multi component minimally 
coupled scalar fields can realize this crossing [14]. 

To study bouncing behavior of the solutions in the Einstein frame with above mentioned 
redefined scalar field, we see that Hubble parameter and scale factor have no possibility to 
run across zero. This behavior has been shown in figure 6 using ansatz presented in the 
previous section. Therefore, in the Einstein frame we have not any possibility to realize 
bouncing solutions with this redefined scalar field. 

4 Summary 

According to existing literature on scalar field dark energy models, a minimally coupled 
scalar field is not a good candidate for dark energy since its equation of state parameter has 
not the capability to realize crossing of the phantom divide line. On the other hand, a scalar 
field non-minimally coupled to gravity in the Jordan frame has the capability to be a suit- 
able candidate for dark energy which provides crossing of the phantom divide line and other 
required facilities. Here we have shown that while a non-minimally coupled scalar field in 
the Jordan frame is a suitable dark energy component with capability to realize phantom di- 
vide line crossing, its conformal transformation in the Einstein frame has not this capability. 
In fact, conformal transformation from Jordan frame to the Einstein frame transforms the 
equation of state parameter of dark energy component to its minimal form with a redefined 
scalar field and in this case it is impossible to realize a phantom phase with possible crossing 
of the phantom divide line. On the other hand, one of the most serious shortcomings of 
the standard cosmology is the problem of initial ( and possibly final) singularity. In recent 
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analysis done within the loop quantum cosmology, the Big Bang singularity is replaced by 
a quintessence Big Bounce with finite energy density of the scalar field. As we have shown 
here, in the Jordan frame with a non-minimally coupled quintessence field, one achieve a 
phenomenologically viable framework to overcome initial singularity with possible realization 
of the bouncing solutions. We have studied the stability of this bouncing solutions too. As 
a result, there are appropriate regions of the uj — uj' phase plane that solutions are stable in 
the Jordan frame. However, the redefined scalar field in the Einstein frame cannot realize 
these bouncing solutions. 
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